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' Abstract 

, The paper proposes, an algorithm to produce novel m-point (for 

any integer m > 2) binary non-stationary subdivision scheme. It has 
been developed using uniform trigonometric B-spline basis functions 
and smoothness is being analyzed using the theory of asymptotically 
equivalence. The results show that the most of well-known binary 
approximating schemes can be considered as the non-stationary coun- 
terpart of the proposed algorithm. 
' Furthermore, the schemes developed by the proposed algorithm has 

, the ability to reproduce or regenerate the conic sections, trigonometric 

polynomials and trigonometric splines as well. Some examples are con- 
, sidered, by choosing an appropriate tension parameter < a < 7r/3, 

' to show the usefulness. 

1 Introduction 

X . 

^ , Subdivision scheme is one of the most important and significant modelling tool to create 



smooth curves from initial control polygon by subdividing them according to some refining 
rules, recursively. These refining rules take the initial polygon to produce a sequence of 
finer polygons converging to a smooth limiting curve. 

In the field of non-stationary subdivision schemes, Beccari et al. [1] presented a 4-point 
binary non-stationary interpolating subdivision scheme, using tension parameter, that was 
capable of producing certain families of conies and cubic polynomials. They also developed 
a 4-point ternary interpolating non-stationary subdivision scheme in the same year that 
generate continuous limit curves showing considerable variation of shapes with a tension 
parameter [2]. A new family of 6-point interpolatory non-stationary subdivision scheme 
was introduced by Conti and Romani [4j. It was presented using cubic exponential B- 
spline symbol generating functions that can reproduce conic sections. Conti and Romani 
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[5] discussed algebraic conditions on non-stationary subdivision symbols for exponential 
polynomial reproduction. 

Since, non-stationary schemes have proven to be efficient iterative algorithms to con- 
struct special classes of curves. One of the important capability is the reproduction or 
regeneration of trigonometric polynomials, trigonometric splines and conic sections, in par- 
ticular circles, ellipses etc. So, in this article, an algorithm has been introduced to produce 
m— point binary approximating non-stationary schemes, for any integer m > 2, using the 
uniform trigonometric B-spline basis function of order m — 1. The proposed algorithm can 
be considered as the non-stationary counterpart of the well known binary approximating 
schemes introduced by Chakin [3] and Siddiqi with his different co-authors [9l [T3l [T^ [T5 t [T6] . 
after setting different values of m in proposed algorithm (for details see table 5.1). More- 
over, the proposed algorithm can also be considered as generalization form of the 2-point 
and 3-point non-stationary schemes presented by Daniel and Shunmugaraj [6]. 

The paper is organized as follows, in section 2 the basic notion and definitions of bi- 
nary subdivision scheme are considered. The algorithm, to produce m— point binary non- 
stationary scheme, is presented in section 3. Some example are considered, to construct the 
masks of 2-point, 3-point and 4-point schemes, in section 4. The convergence and smooth- 
ness of the schemes are being calculated in section 5. Some properties and advantages of 
proposed algorithm are being discussed in section 6. The conclusion is drawn in section 7. 



2 Preliminaries 

In univariate subdivision scheme, following the notion and definitions introduced in [lOj . 
the set of control points {f^ G i? | i G Z} of polygon at k^^ level is mapped to a refined 
polygon to generate the new set of control points {/f "^"^ G i? | i G Z} at the (k + 1)^* level 
by applying the following repeated application of the refinement rule 

= {SaC^)f% = E «S-2,/' V f G Z, ( 2. 1 ) 

with the sequence of finite sets of real coefficients {af'^ , A: > 0} constitute the so-called mask 
of subdivision scheme. If the mask of a scheme are independent of k, namely if for 
all /c > then it is called stationary Sa, otherwise it is called non-stationary S^(k). 

The convergent subdivision scheme, formally denoted by f^~^^ = S^(k)f^, with the cor- 
responding mask af''^ i € Z necessarily satisfies 



j&z jez 

A binary non-stationary subdivision scheme S^(k) is said to be convergent if for every initial 
data G there exits a continuous limit function g G C"^{R) such that 

hm sup|/f -5(2-'=i)| =0, 
and g is not identically zero for some initial data f^. 
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We also recall that two univariate binary schemes S^(k) and are said to be asymp- 
totically equivalent if 

oo 

WSaW - Sfj(k)\\oo < OO, 

fc=l 

where 

1 1 ^aW I loo = max{J2 l4?l> J2 l4i|il} 

Theorem 2.1. The non- stationary scheme S^(k) and stationary scheme Sa are said to 
be asymptotically equivalent schemes, if they have finite masks of the same support. The 
stationary scheme Sa is C"^ and 

oo 

X! ^"^ W^aW ~ 5'a||oo < OO, 
A:=0 

then non- stationary scheme S^ik) is also said to be C™. 

3 The algorithm for approximating schemes 

In this section, an algorithm has been introduced to produce m-point binary non- 
stationary subdivision schemes (for any integer m > 2) which can generate the families of 
(jm-i iijniti]2g curves by choosing a tension parameter < a < vr/S. The algorithm has 
been established using uniform trigonometric B-splines of order n. So, in view of Koch et 
al. [H], trigonometric B-splines can be defined as follow. Let m > n > and < a < ^, 
then Uniform Trigonometric B-splines {T"(x;a)}^^ of order n associated with the knot 
sequence A := {ti = ia : i = 0, 1,2, ...,m + n} with the mesh size a are defined by the 
recurrence relation. 



rd(x;Q) 



1, X € [0, a) 
0, otherwise, 



for 1 < r < n. 



TQ{x;a)= . {sin(x)To ^{x;a) + sin{tr - x)Tq ^{x-a;a)} (3.2) 

sm((r — l)a) 

and TJ{x;a) = Tq{x — ja;a), for j = 1,2, ...,m. The trigonometric B-spline T"(x;a) is 
supported on [tj,tj+n] and it is the interior of its support. Moreover, {TJ^jJL^ are linearly 
independent set on the interval [tn-i,tm+i]- Hence, on this interval, any uniform trigono- 
metric spline S{x) has a unique representation of the form S{x) = J2JLo Pj'^j^ i^'^ o;)' Pj ^ 

To obtain the mask a^(a) = af, i = 0, 1, 2, m — 1 we use the following recreance 
relation, for any value of k, 



af(a) = ro™((m-^-l)^ + ^;^) i = 0, 1, 2, m - 1 (3. 3 
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where yx; ^j, with mesh size (^), is a trigonometric B-sphne basis function of order 
m — 1 and can be calculated form equation (3.2). It can also be observed that the schemes 
produced by the proposed algorithm do not have the convex hull and affine invariance 
properties. Since the sums of the weights of the obtained schemes at kth level are not equal 
to unity. To get sum of the mask equal to unity, the corresponding normalized scheme can 
be obtained (see [6,]). In the following, some examples are considered to produce the masks 
of 2-point, 3-point and 4-point binary approximating schemes after setting m = 2, 3 and 4, 
respectively, in above recurrence relation. 



4 Construction of the schemes 

In this section, some applications are considered to construct the masks of 2-point, 3-point 
and 4-point approximating non-stationary schemes. 

4.1 The 2-point Approximating Scheme 

The linear trigonometric B-spline basis function Tq (^x; , with mesh size (^), can be 
calculated by setting m = 2 in relation (3.3). The 2-point binary non-stationary scheme 
(which is also called corner cutting scheme) with mask af'\a) = af'\i = 0, 1 is defined, for 
any value of k and a g]0, |[, as 



where 



fk+l _ ik) ,k j_ Ak) fk 

hi — % Ji + "^1 Ji+1 

fk+l _ Ak) fk , Ak) fk 

J2i+1 — "l Ji "T "O 



(fc) _ sin 2^ (k) _ sin 



( 4. 4 



sm p sm p 

Theorem 4.1.1. The 2-point binary non- stationary scheme defined above converges and 
has smoothness, for the range a €]0, 
Proof, see [6]. 

Remark 4.1.2. It can be observed that the mask of above 2-point normalized scheme 
converges to the mask of the famous corner cutting scheme introduced by Chaikin [3]. 
Moreover, it can also be observed that the mask of binary 2-point non-stationary scheme 
of Daniel and Shummugaraj [6] can be calculated, after setting m = 2 in (4). Hence the 
proposed scheme can be considered as the generalized form of non-stationary scheme pre- 
sented by Daniel and Shummugaraj and non-stationary counterpart of famous Chaikin's 
scheme [3]. 

4.2 The 3-point Approximating Scheme 

To get quadratic trigonometric B-spline basis function, Tq ^x; with mesh size (^), we 

take m = 3 in equation (3.3). The masks af'\a) = af'\i = 0, 1,2 of the proposed binary 
3-point scheme can be calculated from quadratic trigonometric B-spline function. 
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The 3-point non-stationary scheme is defined, for some value of a €]0, as fohow 

( 4. 5 



J2i — "0 Ji-1 "T "l /« "T "2 Ji+1 



fk+1 _ yK) fk I „V'»'V fk _|_ ^V^V ffc 
J2i+1 — "2 /j-i "T "l /« "T "O 



where 



(k) sm 2^+2 (fe) sm 2k+2 sm 2(1+2 + sm 2k+2 sm 2fc+2 (fc) sm 2*1+2 

'^0 = ^ '■ 2a ^1 ~ ^ a '■ 2a ^2 = ~ a '■ 2a 

sm 2F sm 2F sm 2F sm -^^ sm 2F sm -^^ 

Theorem 4.1.3. The 3-point binary non- stationary scheme defined above converges and 

has smoothness, for the range a e]0, 
Proof, see El- 



Remark 4.1.4. The proposed scheme (4.5) is considered as the generahzed form of the 
non-stationary 3-point scheme developed by Daniel and Shummugaraj in [6j. Furthermore, 
the proposed 3-point scheme (4.5) can be considered as the non-stationary counterpart of 
the stationary scheme |14j . 



4.3 The 4-point approximating scheme 

In this section, a 4-point binary approximating non-stationary subdivision scheme is pre- 
sented and masks of the proposed 4-point binary scheme can be calculated, for any value 
of k, using the relation (4). Where Tq ^x;^^, with mesh size (^), is called the cubic 
trigonometric B-spline basis functions and can be calculated from the recurrence relation 
(3). The proposed scheme is defined, for some value of a s]0, as 



where 



fk+l _ (k) fk j_ „(k) fk I Ak) fk 

J2i — "0 Ji-l "T "1 /j "T "2 /iH 

fk+l _ Jk) fk j_ Jk) fk , (k) fk 

J2i+1 — "3 Ji-i "T "2 /i "T "1 /i 



(k) _ sin^ 2^ 




(4. 6) 



sm 2^ sm 2F sm ^ 



(k) _ sin sin^ + sin sin sin + siiP 
^ sin §T sin || sin || 

(fc) sin 2fc+2 sin 2fc+2 "I" sin 2fe+2 sin gfc+r sin gi+r + sin gfc+r sin 2k+2 



sin ^ sin || sin || 



(A:) _ sin^ 



sm sm 2F sm |jr 



the proposed 4-point scheme can also be considered as the general form of stationary 4- 
point binary approximating scheme, which was introduced by Siddiqi and Ahmad [9]. The 
subdivision rules to refine the control polygon are defined as 



fk+l _ J_fk , 121 fk , 235 fk , 2]_ fk ~) 

J2i — SSi-li-l 384^1 384.'j+l 384-'i+2 I ( A 7 

fk+l _ rj_ k _|_ 235 j-fc I 121 j-fc I _1_ j^fc f I ' 

•'2i+l 384-f^i-i 384-'i 384-'«+l 384-'«+2 J 
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As the weights of the mask of the proposed scheme (4.6) are bounded by the coefficient of 
the mask of the above scheme (4.7). So, we can write as 

(fc) 1 (k) 121 (fc) 235 (k) 27 
° 384' 1 384' 2 384' ^ 384' 

The proofs of a^^'' — > g^^, a^'^^ ^ 02^^^ — )• ||| and 03'^'' ^ can be followed from 
the lemma (5.1.1). 



5 Convergence Analysis 

The theory of asymptotic equivalence is used to investigate the convergence and smoothness 
of the proposed scheme following [8]. Some estimations of 7^^, i = 0, 1, 2, m — 1 are bring 
into play to prove the convergence of the proposed schemes. To establish the estimations 
some inequalities are being considered. 



and 



sma a „ ,7'" 
— r >T for 0<a<b<- 
sm 00 2 

TT 

'csc6 <t CSC t for < 6 < t < — 

2 



sm X TT 

cos X < for < X < — 

~ X 2 



5.1 Convergence Analysis of 4-point Scheme 

To prove the convergence and smoothness of scheme (4.6), estimations of af'\i = 0, 1,2,3 
are being calculated in the following lemmas. 



Lemma 5.1.1. For k > and < a < S 



1 - ^ 1 1 121 ^ (k) . 121 1 



(i) _ < < _ (ii) _ < a «i < 



384 " " - 384cos3(|f) ' ' 384 " ^ - 384cos3(|^) 

..... 235 (fc) 235 1 27 (fc) 27 1 

(ill) < a\ ' < a — (iv) < a\ < a — 

^ ^ 384 -2 - 384cos3(|f) ^ ^ 384 " ^ - 384cos3(|f) 

Proof. To prove the inequality (i) 

(fc) _ sin^(27^) ^ (2^)^ _ J_ 

" sin(f )sin(|^)sin(f ) " (f )(|^)(|^) " 384 



and 



(fe) 



< 



— 23fc+6 



CSC 



—r CSC 

1 



'2a 
.2^, 



CSC 



3a' 
2^^ 
1 



< 



a:^ 1 11 

- 23fc+6 (|f)3cos3(||) - 384^^ 



— 23fc+6 

1 



36 csc^ 



^6a 
.2^ 



6a ^ 
2* ^ 



6 



The proofs of (ii), (iii) and (iv) can be obtained similarly. 

Lemma 5.1.2. For some constants Cq,Ci,C2 and C3 independent of k, we have 

/.X I (fe) 1 1 r-\ I (k) 1 

(1) \-y-sM^<Co^ (") l«i -^1^^123^ 

..... I (fe) 235 1 . I (fc) 27 1 

Proof. To prove the inequality (i) use Lemma (5.1.1), 

384' - 384 cos3a 



< 



< 



1 / 3sin3(||) ^ 
384 1^ 2 cosset 
9a2 1 



64cos3(a) 22fe 

The proofs of (ii), (iii) and (iv) can be obtained similarly. 

Lemma 5.1.3. The laurent polynomial b^i\z) of the scheme {S^cfc)} at the kth level 
can be written as l}i\z) = b''^^\z), where 

6«(z) = 2 { a^\-^ + (af - a^^)z-' + {a^'^ + ^ - a|) + (-^'^^ - af ^ + 4^=) + 4)z 
+ (af ) + af ) - 4^-)).^ + (4'^^ - 4'^)-^ + 

Proof. Since, 

Ak), X {k) -4 , (fc) -3 I (fe) -2 , {k) -1 , (fe) , (fe) , (fe) 2 , (fc) 3 

Therefore using Oq^ + a^^"^ + 02*^^ + 03'^'' = 1, it can be proved. 

Lemma 5.1.4. The laurent polynomial bi{z) of the scheme {S},} at the kth level can 
be written as b\{z) = (^^) where 

b{z) = + 26z-^ + 95 + 140z + 95^^ + 26z^ + z^} 

Proof. To prove that the subdivision scheme {Si^} corresponding to the symbol b{z) is C^, 
we have 



{i + zY 

{z-^ + 22z-' + 1} 



J_|„-2 , r,0_-l 

48 
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Since the norm of the subdivision scheme {Sd} is 

llS'dlloo = max{Y,\d2i\,Y.\^2i+i\} 

r 1 11. 

= max] — , — } < 1 
So in view of Dyn [lOj, the stationary scheme is C^. 

Theorem 5.1.5. The 4-point non- stationary scheme defined in Eq. (4-6) converges 
and has smoothness, for the range a €]0, 

Proof. To prove the proposed scheme to be C^, it is sufficient to show that the scheme 
corresponding to the symbol b^{z) is (in view of Theorem (8) given by Dyn and Levin 
[8]. Since {Sf,} is by Lemma (5.1.4). So, it is sufficient to show for the convergence of 
binary non-stationary scheme (4.6) that, 

oo 

J22^''\\Sb(^)-Sb\\oo<00 
k=0 

Where, 



Sfjik) - Si 



felloo = max 



+2j ~ 



0,1} 



Following Lemmas (5.1.3) and (5.1.4), we have 



(fc) 
2i 



52i 



(fc) 


1 





~ 384 


(fc) 


1 





~ 384 



and similarly, it may be noted that 



ik) 



384 



+ 2 



ik) 



+ 4 
+ 4 

121 

~ 384 



(fc) (fc) (fc) 95 



a\ + a 



384 



Sk) 



+ 2 



121 
384 



ik) 



+ 4 



ik) 



235 
384 



+ 6 



27 
384 



ik) 



27 
384 



From (i), (ii), (iii) and (iv) of lemma (5.1.2), ^ 23^=14'^ -3I4I < 00, „ 23^=1 af^-lgl < 



235 
384 



< 00 and Yl,T=o '^^^\o!''^^ — ^| < 00. Hence, it can be written as 



k=0 

Thus by the Theorem 2.1, c'^^^z) is as the associated scheme c{z) is . Hence, 
proposed scheme is C^. 
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6 Properties and advantages of algorithm 



In this section, some properties hke unit circle reproduction property, symmetry of basis 
function and some other advantages of the proposed algorithm are being considered. 

6.1 Reproduction of unit circle 

It can be observed that certain functions like cos(q!.) and sin(a.) can be reproduced 
by the proposed non-stationary schemes. In particular, if a set of equidistant point f-' = 
cos((z — i)a) and a = then the limit curve is the unit circle (see also figure 2). Similarly, 
sm{{i — \)a) reproduces sin(a.). 

Proposition 6.1. The limit curves of the scheme (4.4) reproduces the functions cos(a.) 
and sin(Q;.) for the data points = cos((i — 5)^) and ff^ = sin((i — ^)-^), respectively. 
In other wards we have to show that for > 

(i) ff^ = cos((i — 5)^), we have for < i < 2^n as 

= cos (2z^) and /.^^^ = cos ((2z + 1)^) 

(ii) Similarly, for //^ = sin((i — ^)^), we have for < i < 2''n as 

fi+' = sin (2i^^ and f^^' = sin ((2z + 1)^^) 

Proof: For any initial data of the form = cos((i — ^)a), it can be followed 

fl _ „{0)fO , (o).o 

sm(f) 1 .^sm(f) ,,.^1. , ,. . 

= ^ cos u a) H ^ cos z H — )a) = cos(zQ;) 

sina 2' ' sina 2^ ^ ^ ' 

Similarly, it can be followed for fe*'* refinement 

a 



/I, = cos {{21) 



Analogously, it can also be proved. 



f2i+l = COS (2i + 1 



' 2^'+i 



Proof of part (ii) is similar. So, the scheme (4.4) can reproduce sin(a.) for the initial data 
/0 = sin((i-i)a). 

Proposition 6.2. The limit curves of the scheme (4.5) reproduces the functions cos(a.) 
and sin(a.) for the data points = cos((i + ^)^) and //^ = sin((i + ^)-^), respectively. 
It can be proceed on same way. If a set of equidistant point = cos((i — ^)a) and 
a = ^ can be chosen by the scheme (4.5), then the limit curve is the unit circle. Similarly, 
sin((z — ^)a) reproduces sin(a.). 
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6.2 Symmetry of Basis Limit Function 

The basis limit function of the scheme is the limit function for the data 



In order to prove that the basis limit function is symmetric about the Y-axis. 



1, i = 0, 
0, i^O 



Theorem 6.1 The basis limit function F is symmetric about the Y-axis. 

Proof. The symmetry of basis limit function can be followed on the same pattern following 



6.3 Special cases 

It can be observed that the binary subdivision schemes presented in [3l O [U [13l [HI [151 US] 
are either the special cases or can be considered the non-stationary counterpart of the 
stationary schemes. (See also Table 1 and Table 2). 

• The hmiting curves can be obtained after taking m = 2 in proposed algorithm. 
The obtained curves of scheme (5), taking a = coincide with the limit curves of 
the famous corner cutting scheme of Chaikin |3J. 

• After setting m = 2 and m = 3 in proposed algorithm, the mask of 2-point and 3-point 
approximating non-stationary schemes, developed by Daniel and Shunmugaraj [6j, can 
be obtained. 

• The limit curves of 3-point stationary approximating scheme, introduced by Siddiqi 
and Ahamd HI], coincide with the limit curves of proposed 3-point scheme, after 
setting m = 3 and a = 

• The limit curves of proposed 4-point scheme coincide with the limit curves obtained 
by the scheme [9j , for m = 4 and a = . 

• The limit curves of 5-point stationary approximating scheme introduced in |16j matched 
with the limit curves of proposed scheme, for setting m = 5 and a = 

• The limit curves of 6-point stationary approximating scheme introduced in matched 
with the limit curves of proposed scheme, for setting m = 6 and a = 

• The proposed algorithm can also be considered as the non-stationary counterpart of 
m-point scheme developed by Siddiqi and Younis [15j . 
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Table 1: The proposed algorithm can be considered as the non-stationary counterpart of 
the following existing binary approximating schemes, for different values of m. 



Setting m 


Scheme 


Type 


Continuity 


Counterpart 


2 


2-point 


Stationary 




scheme ^ 


3 


3-point 


Stationary 




scheme [14] 


4 


4-point 


Stationary 




scheme [S] 


5 


5-point 


Stationary 




scheme [TB] 


6 


6-point 


Stationary 




scheme [13] 


m 


m-point 


Stationary 


1 


scheme [TS] 



Table 2: The following existing binary non-stationary approximating schemes can be con- 
sidered the special case of proposed algorithm, for different values of ra. 



Setting m 


Scheme 


Type 


Continuity 


Coincides with 


2 


2-point 


Non-stationary 




scheme 


3 


3-point 


Non-stationary 




scheme [B] 



7 Conclusion 

An algorithm of m— point binary approximating non-stationary subdivision scheme 
(for any integer m > 2) has been developed which generates the family of (7™"^ limiting 
curve, for < a < vr/S. The construction of the algorithm is associated with trigonometric 
B-spline basis function. It is also evident from the examples that the limit curves of the 
proposed schemes coincide with the schemes presented in [S] [6l [HI [9l [131 US]- So, the 
proposed algorithm, for different values of m, can be considered as the generalized form 
of the scheme [6] and non-stationary counterpart of the stationary schemes [3l [9l [131 El 
I16j . Moreover, the schemes produced by the algorithm can reproduce or regenerate the 
trigonometric polynomials, trigonometric splines and conic sections as well. 
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Figure 1: The continuous lines show the behavior of Umit curves (after three subdivision 
steps) of the proposed scheme in blue, green, black and red colours, respectively, after 
setting m = 2,3,4 and 5. The blue, green, black and red coloured limit curves coincide 
with the schemes presented in [3], [6], [9] and [16], respectively. The dotted lines represent 
the control polygon 
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Figure 2: The green, blue and magenta coloured continuous lines represent the limit curves 
of scheme [3], [H] and [9], respectively in (i). The blue coloured continuous line represents 
the limit curve of proposed scheme in (ii).The reproduction of the unit circle and ellipse 
(the limit curves in continuous lines are obtained after three subdivision steps). The red 
coloured broken lines represent unit circles and ellipses and black coloured dotted lines 
represent control polygon. 
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